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Abstract:  Utilization  of  controllable  ferroelectric  and  ferromagnetic  layers  coat¬ 
ing  a  conducting  object  to  provide  an  attenuation  capability  against  electro¬ 
magnetic  interrogation  is  discussed.  The  problem  is  formulated  as  a  differential 
game  and/or  a  robust  optimization.  The  scattered  field  clue  to  interrogation 
can  be  attenuated  with  the  assumption  of  an  uncertainty  in  the  interrogation 
wave  numbers.  The  controllable  layer  composed  of  ferromagnetic  and  ferroelec¬ 
tric  materials  [9,  10]  is  incorporated  in  a  mathematical  formulation  based  on 
the  time-harmonic  Maxwell  equation.  Fresnel’s  law  for  the  reflectance  index  is 
extended  to  the  electromagnetic  propagation  in  anisotropic  composite  layers  of 
ferromagnetic  and  electronic  devices  and  used  to  demonstrate  feasibility  of  con¬ 
trol  of  reflections.  Our  methodology  is  also  tested  for  a  non-planar  geometry  of 
the  conducting  object  (an  NACA  airfoil)  in  which  we  report  our  findings  in  the 
form  of  reduced  radar  cross  sections  (RCS). 


1  Introduction 


In  this  paper  we  discuss  an  optimal  attenuation  problem,  i.e.,  we  attempt  to  maximize 
attenuation  capabilities  of  interrogating  signals  by  utilizing  a  controllable  dielectric  layer 
on  the  surface  of  a  conducting  object.  The  objective  of  the  interrogator  is  to  detect  and 
identify  the  location  and  shape  of  the  conducting  object  based  on  the  scattered  field  from  an 
interrogation  incident  field,  i.e.,  the  solution  of  an  inverse  scattering  problems  [2,  6].  In  the 
plane  wave  case  the  incident  electromagnetic  (EM)  field  has  the  form  (E^\  H^)  elk'x  and  the 
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interrogator  has  control  over  the  wave  numbers  k.  The  attenuation  problem  is  to  minimize  or 
diminish  the  detection  capability  of  interrogation  by  either  grating  [13]  or  fabrications  on  the 
object  surface.  In  this  paper  we  consider  the  utilization  of  thin  controllable  dielectric  surface 
layers  on  the  object  as  a  method  for  achieving  the  attenuation  capability.  Here  “controllable” 
means  that  we  have  the  capability  of  adjusting  the  material  properties  of  the  surface  layers 
parametrically  and/or  dynamically.  The  technical  ideas  developed  here  also  have  potential 
to  aid  in  the  design  of  medical  shields  employed  to  protect  parts  of  an  irradiated  target  or 
to  focus  radiation  to  pinpoint  specific  regions  of  the  target. 

In  these  initial  investigations,  we  investigate  an  “active”  design  case  in  which  one  determines 
values  of  the  dielectric  permittivity  and  magnetic  permeability  of  the  controllable  layer  in 
order  to  attenuate  reflections.  From  a  control  theoretic  viewpoint  this  is  a  “passive”  or 
open  loop  control  strategy.  But  our  efforts  here  lay  the  foundations  for  “active”  or  closed 
loop  control  strategies  in  which  one  combines  controllable  layer  dynamics  with  a  sensor  for 
incoming  interrogating  signals  to  develop  real  time  feedback  controls  for  adaptive  choice  of 
the  permittivity  and  permeability  of  the  controllable  layer. 

To  investigate  feasibility,  we  first  formulate  the  problem  as  a  differential  game.  For  example, 
we  assume  the  time-harmonic  incident  EM  plane  wave  is  impinging  on  the  surface  at  z  =  0 
and  we  control  the  effective  dielectric  constant  e  of  the  surface  layer  on  top  of  the  conducting 
material.  The  scattered  field  due  to  the  interrogation  can  be  evaluated  based  on  the  time- 
harmonic  Maxwell  equations  [6].  In  the  case  when  the  dielectric  constant  is  homogeneous  in 
the  horizontal  directions  (planar  geometry)  the  reflectance  index  R  =  R{k ,  e)  is  determined 
by  Fresnel’s  law  (see  [14]  and  Section  4).  Thus,  the  problem  of  nullifying  the  scattered  held 
can  be  cast  as  the  minimization  of  the  scattering  EM  wave  in  terms  of  |i?|2,  i.e., 

min  max  |i?|2  (1.1) 

ee(2  iieK 

where  Q  is  a  set  of  admissible  dielectric  constants  and  K  is  a  set  of  possible  interrogation 
wave  numbers.  In  order  to  determine  the  admissible  set  Q  we  must  describe  controllable 
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mechanisms  of  the  dielectric  layer.  Thus,  the  problem  of  nullifying  the  scattered  held  can 
be  formulated  as  min-max  problem  of  minimizing  the  largest  reflectance  by  interrogations 
over  k  G  K,  over  all  possible  designs  (e  G  Q  in  this  case).  This  game  theoretic  formulation  is 
used  in  many  other  design  problems.  It  does  not  assume  any  information  on  the  uncertainty 
of  interrogations  and  thus  it  may  lead  to  a  conservative  design.  An  alternative  formulation 
can  be  given  in  a  more  robust  form,  i.e., 

min  [  \R\2  dJC(k)  (1.2) 

£eQ  JK 

where  fC  is  a  probability  distribution  function  on  the  wave  numbers  k.  This  formulation 
then  needs  information  about  the  distribution  function  /C  of  the  interrogating  plane  wave. 
As  demonstrated  in  Section  2,  better  knowledge  of  the  distribution  function  greatly  improves 
the  performance  of  minimizing  the  scattered  held. 

After  demonstrating  attenuation  capabilities,  we  turn  to  the  general  case  of  a  non-planar 
conducting  medium  with  a  controllable  coating  layer.  The  far  held  pattern  F{6)  of  the 
scattered  held  (see,  [6]  and  Section  5)  is  then  a  function  of  the  wave  number  k  of  the 
incident  plane  wave  and  the  material  properties  (e,/i)  of  the  controllable  layer,  i.e., 

F{0)  =  U{k,{e,n)\Q),  0  <  #  <  27t.  (1.3) 

One  can  select  the  performance  index  J(k,(e,n))  =  Q(F)  to  perform  specihc  alternations 
of  the  scattered  held,  which  of  course  depend  on  the  inverse  techniques  employed  by  the 
interrogator.  Here  Q(F)  is  some  performance  index  for  the  far-held  pattern  F{6).  We 
investigate  optimal  radar  cross  sections  (RCS)  for  one  class  of  such  problems. 

We  note  that  the  existence  of  solutions  to  a  general  min-max  problem  is  guaranteed  under  the 
condition  that  given  e  G  Q  the  value  function  V(e)  =  sup  J(k,e )  is  lower  semi-continuous, 

keK 

which  is  typically  satisfied  under  very  mild  conditions  (e.g,  see  [12])  when  Q  is  compact. 
The  saddle  point  property  of  a  solution  pair  (fc0,eo): 

J(k,e0)  <  J(k0,e0 )  <  J(k0,e )  for  all  e  G  Q,  k  G  K  (1.4) 
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holds  locally  if  the  Hessian  of  J(k,e )  is  hyperbolic  at  ( fcoAo ).  The  existence  of  solutions  to 
the  robust  formulation  then  simply  follows  from  the  continuity  of  J(k,e o)  with  respect  to 
(k,  e)eKx  Q. 

A  brief  outline  of  our  presentation  here  is  as  follows.  In  Section  2  we  consider  the  exact 
problem  for  the  planar  geometry  and  present  example  numerical  calculations  for  optimal 
design  of  the  dielectric  layer  based  on  a  robust  formulation  and  demonstrate  the  feasibility 
of  this  approach  and  the  effectiveness  of  the  design.  In  Section  3  we  discuss  a  controllable 
layer  composed  of  ferromagnetic  and  ferroelectric  materials  as  proposed  by  How  and  Vittoria 
[9,  10].  This  composite  model  is  designed  so  that  a  control  mechanism  for  the  material 
properties  of  the  layers  can  be  achieved  in  both  a  parametric  and  a  dynamic  manner.  In 
Section  4  we  mathematically  formulate  the  forward  problem  for  the  controllable  composite 
layers,  including  the  tensor  permeability  in  the  ferrite  layer,  by  calculating  the  scattered  field 
R  as  a  function  of  the  interrogating  wave  and  the  near  surface  composition.  We  consider  the 
time- harmonic  case  with  plane  wave  interrogations  of  the  planar  geometry  (i.e.,  the  composite 
layers  are  homogeneous  in  ( x:y )).  In  this  case  we  construct  the  plane  wave  solution  and  an 
analytic  expression  for  the  reflectance  index  R.  Another  important  feature  of  our  formulation 
is  the  possible  identification  of  the  interrogating  wave  in  terms  of  its  distribution  /C.  Since 
the  plane  wave  calculation  also  yields  surface  currents  as  an  explicit  function  of  the  incident 
interrogations,  the  surface  current  measurements  can,  in  principle,  be  used  to  identify  the 
distribution  /C  of  the  interrogations  k.  In  Sections  5-7  we  present  results  for  our  formulation 
when  applied  to  a  non-planar  geometry  by  considering  the  NACA0012  airfoil  [15].  In  general 
(and  in  particular  in  this  case)  we  do  not  have  an  analytic  expression  for  the  far  field  pattern 
F  and  thus  we  use  a  numerical  computation  of  the  scattered  field.  In  this  case,  our  numerical 
computations  for  the  scattered  field  F  are  implemented  using  the  finite-element  method  (in 
Section  6).  Our  numerical  findings  are  presented  for  an  optimal  homogeneous  coating  layer. 
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2  Feasibility  study 


In  this  section  we  first  demonstrate  the  feasibility  of  our  approach.  We  consider  an  incident 
parallel-polarized  (TEX)  plane  wave  H  =  (H^\  0,  0)elk  x  impinging  on  the  interface  of  the 
first  and  second  layers  at  z  =  0,  as  depicted  in  Figure  1.  The  interface  between  the  second 
and  third  layers  is  located  at  z  =  —  d,  where  the  third  layer  is  a  perfect  conductor.  We  control 
the  effective  dielectric  constant  e®  of  the  second  layer  coating  the  conducting  material.  By 
Fresnel’s  law  (see  [3,  14]  and  the  discussion  in  Section  4)  we  have  that  the  reflection  coefficient 
or  reflectance  index  is  given  by 


R 


=(2)fc*1}  ~  e(1)#  +  e_2  ik^d 


d2)/^  +  e<4)/c(2') 


=(2)fc*1}  -  e(1)^2)  -2 ik[2)d 


1  + 


-P'>ki1'>  +  ed)/cl2') 


(2.1) 


where 

=  y  \j 1  -  sin2  <Po 

kf  =  ^  sin2  ip0. 

A  * 


(2.2) 


(i) 


Here  tpo  is  the  incident  angle  with  respect  to  the  normal  to  the  surface  (tan</?0  =  ~7j\,  k^  = 


rvz 

0)  and  A  is  the  wavelength  of  the  incident  wave.  We  note  from  (2.1)  that  R  depends  on  the 
ratio  /e^  and  hence  without  loss  of  generality  we  may  normalize  the  parameters  so  that 
e(!)  _  i  The  reflectance  index  R  is  a  function  of  (A,<^o),  the  normalized  dielectric  constant 
e^2\  and  the  thickness  d  of  the  surface  layer.  We  assume  that  d  is  positive  and  fixed.  We 
parameterize  the  incident  wave  in  terms  of  (A,</?o)- 


In  Figures  2  and  3  we  depict  the  robustness  of  the  optimal  solution  by  plotting  the  reflectance 
intensity  \R\  as  a  function  of  the  incident  angle  and  the  normalized  thickness/wave  length 
ratio,  defined  as  a  —  In  Figure  2  we  assume  that  the  uncertainty  in  wave  numbers  is  only 
due  to  uncertainty  in  the  incident  angle  ipo,  which  is  uniformly  distributed  on  the  interval 
[36,  54]  degrees,  and  graph  the  intensity  \R\  corresponding  to  the  optimal  dielectric  constant 
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Figure  1:  Schematic  representation  of  the  reflection  of  a  plane  wave  incident  with  angle  <j>o  on  a  planar 
three  layer  stack.  The  top  two  layers  are  dielectric  media.  The  third  layer  is  a  perfect  conductor. 

e*-2'  =  1.008  +  .693?  as  a  function  of  the  incident  angle.  The  integration  of  \R\ 2  over  ipo  is 
performed  using  Simpson’s  rule.  The  intensity  of  the  reflection  is  well  attenuated  over  the 
uncertainty  interval  [36,  54]  degrees. 


Optimal  1 . 008+0. 693i  for  a=0.5  and  on  [36  54]  degree 


Figure  2:  The  reflected  intensity  \R\  as  a  function  of  the  incident  angle 
ifQ.  The  uncertainty  interval  is  [36,  54]  degrees. 

Next,  we  assume  that  there  is  uncertainty  in  both  the  incident  angle  ip 0  and  the  normalized 
thickness/ wave  length  ratio  a,  which  are  uniformly  distributed  on  a  rectangle  [36,54]  x 
[0.3, 0.7].  In  Figure  3  we  plot  the  intensity  \R\  corresponding  to  the  optimal  dielectric 
constant  e ^  =  1.4309  +  1.0724?  for  several  sampled  incident  points  in  the  frequency.  A 
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reasonable  attenuation  over  the  uncertainty  box  is  obtained  in  this  example.  It  is  clear 
that  the  performance  depends  on  the  quality  of  the  information  on  the  distribution  of  the 
interrogating  wave,  as  demonstrated  by  the  better  attenuation  results  in  the  first  case  (Fig. 
2)  than  in  the  second  case  (Fig.  3). 


Optimal  1 .4309+1 .0724i  for  a  in  [0.3  0.7]  and  on  [36  54]  degree 


Figure  3:  The  reflected  intensity  \R\  as  a  function  of  the  incident  angle 
ipo  and  the  normalized  thickness /wave  length  ratio  a.  The  uncertainty 
box  is  [36,  54]  x  [0.3,  0.7] 
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3  Controllable  sublayers  composed  of  ferromagnetic  and 
ferroelectric  materials 


In  this  section  we  describe  an  experimental  device  that  can  be  used  to  control  the  dielectric 
permittivity  (and  magnetic  permeability)  in  a  coating  layer  as  discussed  in  the  previous 
sections.  In  Figure  4  we  present  a  schematic  of  the  configuration  of  an  active  reflecting  device 
proposed  and  investigated  experimentally  by  How  and  Vittoria  in  [9,  10].  The  reflector 
contains  a  ferrite  layer  and  a  ferroelectric  layer  as  constituents.  The  permanent  magnet 
provides  a  common  magnetic  bias  so  that  the  ferromagnetic  resonance  (FMR)  condition  can 
be  readily  achieved  and  thereby  facilitate  sensitive  magnetic  tuning  by  the  local  Helmholtz 
coils.  The  dielectric  properties  of  the  ferroelectric  layer  are  controlled  through  the  ground 
plane  bias  held.  The  purpose  of  this  rehector  design  is  to  provide  phase  and  impedance 
control  of  the  composite  layers  so  that  nullification  and  alteration  of  the  scattered  wave  can 
be  achieved  in  the  response  to  an  incident  interrogating  EM  wave.  The  integrated  circuits 
are  designed  so  that  the  tuning  sensitivity  of  the  device  is  enhanced.  The  key  element  of  the 
device  is  that  the  material  properties  y(H)  and  e(E)  of  the  composite  layers  are  controllable 
in  terms  of  the  magnetic  mean  in  the  ferrite  layer  and  the  electric  mean  in  the  ferroelectric 
layer,  and  thus  can  support  agile  frequency  attenuation. 


The  most  important  device  characteristic  of  the  ferrite  in  our  study  [9,  10]  is  that  the 
magnetic  permeability  y  is  a  tensor,  so  that,  due  to  the  gyromagnetic  effect,  EM  propagation 
in  the  ferrite  is  anisotropic  in  the  presence  of  a  dc-bias  magnetic  field  [11,  19].  For  a  ferrite 
magnetized  in  the  y  direction  with  damping  and  no  demagnetization,  the  permeability  tensor 
is  given  by  [16] 


h  = 


y  0  —in 
0  y0  0 


\  OJOJ 

h  =  ho(l+^2 - K  = 


LVq  —  U)2' 


where  ~u  ~  U  (3.1) 

Ujm  =  47T7  Mz,  U0  =  (V0  +  if  t  ,  cu0  =  7  H0, 


in  0  /i 

u>0  is  the  precession  frequency,  H0  is  the  impressed  d.c.  magnetic  held,  7  is  the  gyromagnetic 


ratio,  Mz  is  the  saturation  magnetization,  and  r  is  the  relaxation  time. 


The  ferrite  device  is  most  useful  if  it  operates  near  the  FMR  frequency  Uq  so  that  the  rapid 
change  in  magnetic  permeability  can  be  effectively  utilized,  either  to  obtain  frequency-tuning 
capability  or  to  remove  the  degeneracy  between  modes  [9,  10]. 
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Figure  4:  Composite  sublayers  comprising  an  active  reflector  device. 


4  Plane  wave  solution 


We  next  discuss  a  plane  wave  solution  as  it  interacts  with  a  ferrite  layer.  Due  to  the  tensor 
magnetic  permeability  the  electric  and  magnetic  modes  are  coupled  in  the  ferrite  layer.  In 
this  section  we  present  the  detailed  calculations  for  constructing  the  fundamental  solution 
in  the  ferrite  layer.  A  similar  calculation  can  be  carried  out  for  the  ferroelectric  layer. 
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The  time-harmonic  Maxwell’s  equations  are  written  as 


V  x  H  =  iueE 

V  x  E  =  —iujlH 
V  •  (p,H)  =  0 

V-E  =  0. 


For  the  ferrite  layer  the  permeability  [X  is  the  tensor  defined  by  (3.1),  and  thus  the  first  three 
equations  in  (4.1)  can  be  written  as 

(  _  Air  \ 

%  if  Ex\ 

—Hx  -  —Hz  =  iue  [  Ey  ,  (4.2) 

\EJ 

\  dx  y  dy  x  ) 


dr,  d  r,  \ 

ay  dz 

d^  d  ^ 

—Ex  -  —Ez  =  -iu 
dz  dx 

^-Ey  -  ^-Er 
dx  y  dy  ) 


yHx  —  inE[z 
yoHy 

zkHx  -(-  yHzl 


-,l  ^dxHx  +  dzHz)  ~  in{dzHx  dxHz)  +  ^  dyHy' 


Taking  the  cross-product  of  (4.2)  and  using  (4.3),  we  obtain  (^-component) 


.  if  d2  TT  d2  u  d  f  d  TT  d  TT 

xujfiQ  Hy  .  too  o  o  **  y  q  (  ^  ** x  ^  ** z 

Hue  \oxz  ozz  ay  \ox  oz 


Using  the  ^/-component  of  (4.2), 


iueEr,  = 


d_H 

r\  -H-X  r\  -H-Z1 

OZ  OX 


in  (4.4)  and  substituting  into  (4.5),  we  have 


92  „  ,  Mo  9 2 


cjKe  d 


dx*Hy+dz*Hy+  »  dy*Hy+U  ^eHy~  //  dyEy' 
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The  components  of  (4.3)  can  be  manipulated  to  yield  the  following  equations 

d  d  d  d 

n  (— Ez  -  — Ey )  +  m  (— Ey  -  —  Ex)  =  -  k2)Hx 

dy  dz  dx  dy 

d  7 _ >  d 

~dz^x  ~dx^z  Hy 


d  c)  d  3 

-in  (— £*  -  —Ey )  +  /i  (— Ey  -  — -E1*)  =  -iixi/j2  -  n2)  Hz. 


'  dy  dz 


dx  dy 


Using  the  ^/-component  of  (4.2)  and  the  first  and  last  equations  of  (4.8),  one  can  obtain 

d  (  ,  d  ^  d  ^  x  <9 

Wz  ^dy  z  ~  d~z  y)  +  lK  {dxEy  ~  dy  x)) 


dx  (  lK^d  yEz  dzEy)  +  *  ^dxEy  dyEx^)  ^ 2  Ey 


This  can  be  rearranged  to  obtain 

( d2  d2  \ 
_/i  (iu2  y  + 


+o>2/  m(a^*  dxEz)S)  ~uj2e^2  K^Ey 

Using  this  result  with  the  ^/-component  of  (4.8)  and  the  fourth  equation  of  (4.1),  we  have 

d2  d2  d2  u2e  (/i2  —  k2)  lu/j>ok  d 

fa2  y  +  dfEy  +  d^y+  Ji  y  =  ~iT  dyHy '  ^  “ 


From  (4.2) 


If  we  define 


—  U  -  —  H 

dy z  dz  y 

ax  a?/ 


£/  —  Ez  =t  i  Ex,  iiz  =  Hz  it  i  Hx, 


then  it  follows  that 


-i\7+Hy  +  i— 


(4.10) 
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where 


V±0=—  <j>±i  —  c 
dz  ox 


From  (4.3) 


and  thus 


H  Hx  —  %kHz 
iiHz  +  in  Hx 


(n  +  k)  H 4 
(p  —  k)  H 


d  d 
dyEz  ~  dzEy 

^-Ey  -  ^~EX 
dx  y  dy 


-iV+Ey  +  i-^-E+ 
dy 


iV  Ev  —  i—E 
dy 


From  (4.10)-(4.11) 


d2  d 

-cu2e  (fj,  ±k)H±-  ±  cue  V± Ey  +  —  =  0, 

dy2  dy 


and  we  have 


d2  d 

-c A  (n  ±k)E±-  —E±  t  c o(ii  ±  k)  V± Hy  +  — ' V±. Ey  =  0. 

dy 2  dy 

Let  E  and  E[  be  the  (partial)  Fourier  transform  (in  (. x,y ))  of  E  and  H,  i.e., 


(4,11) 


E(kx,  hy ,  z') 

Z7T 


E(x,  y,  z)e 


—ikx  x—iky  y 


dxdy. 


Then  (4.9)  and  (4.7)  can  be  written,  respectively  as 


d2  ~  II2  —  K2  „  K  ~ 

E,  -  (kl  +  kl-  )  E,  =  -iu^ky- H, 


q  2  Ey  ( hx  T  hy  cu^/xq^)  hdy  xl uahy  '  Ey 


and  (4.10)-(4.11)  as 


—  iky  V1*1  Ey  ±  cu(/x  ±  k)  S7±Hy 


E±  = 


H±  = 


k2  —  cu2e(p  ±  k) 

-iky  \/±Hy  =F  cue  V±Ey 
h2  —  cu2e(p  ±  k) 


(4.12) 


(4.13) 
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where 


V±  =  ^  T  kx. 

oz 


Thus  the  electric  and  magnetic  modes  are  coupled  by  (4.12).  Next  we  find  the  fundamental 
solution  within  the  ferrite  layer.  We  look  for  a  solution  satisfying 


Ey  —  i  f  i  Hy. 


(4.14) 


Then  rj  must  satisfy 


2  2  h2  -  «2  ,  «  ho ,  2  2  ,  « 

k  —  u  e - utiQky —  =  — k  —  uj  fi 0e  —  uekv—rj. 

y  H  nrj  /i  y  n 


(4.15) 


Equation  (4.15)  has  two  roots,  which  we  designate  as  rj+  and  r/_ .  Then  Hy  must  satisfy 

~  (kl  +  -  u2ii0e  -  r/±  ueky^)  Hy  =  0. 

Thus,  Ey  and  Hv  can  be  written  in  the  form 


Ev  =  irj+  (A1  eik+z  +  B1  e~ik+z )  +  ir)_  (A2  eik~z  +  B2  e~ik~z ) 
Hv  =  A!  eik+z  +  B1  e~ik+z  +  A2  eik~z  +  B2  e~ik~z, 


(4.16) 


k±  =  Ju)2/i0e  +  ri±ueky—  —  k2  —  —k2. 


Formulae  (4.12)-(4.16)  are  given  in  [11,  19]  without  detailed  derivations.  These  results  allow 
us  to  construct  the  plane  wave  solution  in  a  ferrite  layer. 

Moreover,  in  a  general  dielectric  (including  ambient)  medium  the  fundamental  solutions  (in 
the  partial  Fourier  domain  formulation)  are  given  by 


(Ex,  0, 0)eikzZ,  (0 ,-—Ex,  ^Ex)eikzZ 

U)[l  CO/J, 

for  a  perpendicular  polarized  (TM_r  mode)  incident  wave  and 


(4.17) 


(0,  hfl -^Hx)e‘k“,  (Hx, 0,0)e“V 

cue  cue 


(4.18) 
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for  a  parallel  polarized  (TES  mode)  incident  wave,  where  Ex  and  Hx  are  constants,  when 
kx  —  0  and  ky  are  fixed. 

Now  we  consider  the  case  when  a  ferrite  layer  with  thickness  d  on  a  perfectly  conducting  is 
impinged  upon  by  the  parallel  polarized  incident  wave  (4.17)  H^\x,y,  z)  =  (H^\  0,  0)elk  x 
(its  partial  Fourier  transform  is  (H^\  0,  0)eikzZ  )  with  kx  =  0  and  ky  fixed.  The  transmitted 
wave  in  the  ferrite  layer  is  generally  not  TEx-mode  alone  and  thus  has  a  nontrivial  Ex\ 
Hence,  for  a  given  incident  wave,  the  reflected  wave  (E^r\  H^)e~‘lkzZ  in  the  ambient  layer 
is  a  linear  combination  of  the  two  fundamental  solutions  of  the  form  (4.17)-(4.18)  (with 
kz  =  —  kz)  and  the  constant  weights  (E^\  H^)  for  the  reflected  wave  can  be  determined  by 
the  system  of  equations 

'  ffW  +  fff  =  ff»(0) 

—  (o  +  4r>)  =  ffW(o) 

U/J,  y 

<  0  +  E^  =  E®(0)  (4.19) 

-  («?  -  ffb)  =  Ef( 0) 

cue  y 

,  Ei‘\~d)  =  EV(-d)  =  0, 

for  (H1?  Bi,  A2,  E>2,  E^\  H^).  The  hrst  two  equations  impose  the  continuity  of  H  compo¬ 
nents,  the  next  two  equations  impose  the  continuity  of  E  components  at  z  =  0  (the  interface 
between  the  ambient  and  the  ferrite  layer),  and  the  last  equation  enforces  the  perfectly  con¬ 
ducting  boundary  conditions  at  z  =  —d.  Moreover,  the  induced  surface  current  J  due  to  the 
incident  wave  is  given  by 

J  =  (. Hy(-d ),  ~Hx{-d),  0)  =  n  x  H. 

This  construction  procedure  can  be  readily  extended  to  the  case  of  the  composite  of  sublayers. 
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4.1  Dielectric  case  (Fresnel’s  Law) 


In  this  section  we  consider  the  dielectric  layer  (fj,  —  /i0,  n  =  0  and  thus  rj  =  0)  with 
thickness  d  on  the  perfectly  conducting  medium  and  show  that  equation  (4.19)  reduces  to 
usual  Fresnel’s  law  for  the  parallel  polarized  (TE_r)  incident  wave.  In  this  case  we  have 
E M  =  E®  =  =  H®  =  0  and 

kf]  =  ^/cu2^1)  -  k2 

k+  =  k-  =  kf^  =  0e(2)  —  k2. 


In  the  dielectric  layer  we  have  (in  the  partial  Fourier  domain) 

H®  =  Hf  eik*)z  +  H®  e~ik™z 
„(2) 


k\ 


Ed)  =  <H 

V  --(2)  V 


r(0  _  jjd)  e-*fc)2) 


cue 


where  H±*  are  constants.  Thus,  (4.19)  becomes 


kb 


Hf  +  H?  =  Hf  +  (Hf  -  Hi'>)  = 


1.(2)  7,(2) 

kz  R(t)  _  kz_  u{t) 


cue 


e-ikPdHd)  _  eik^dHd)  =  0 


cue 


(2) 


cue 


(2) 


H' 


Hence  we  obtain  Fresnel’s  law 


Hir)  = 


A2)kW  _  €(l)fc(2) 

d2)/cl1')  +  ed)k^ 


+  e 


-2  iki2)d 


t{2)k{z]  ~  e(1)^2)  -2ikwd 


H®. 


1  + 


e(2)/ci1')  +  ed)/ci2-) 


4.2  Numerical  tests 

In  this  section  we  demonstrate  the  feasibility  of  using  the  controllable  property  of  the  ferrite 
layer  to  attenuate  reflections.  We  select  e0  =  1,  /i0  =  10  and  d  =  .5  (normalized).  The 
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magnetic  permeability  fx  of  the  ferrite  layer  is  parameterized  by 

H  =  H o(°-3+^)>  «  =  /^o(-0.01  +  ^),  for  1  <  /?  <  100 

where  a  parameter  (3  plays  a  role  of  tuning  the  frequency  um  in  (3.1).  In  Figure  5  we  depict 
\{E^\H^)\  as  a  function  of  (3  for  the  three  different  incident  angles  0O  =  40°,  45°,  50°. 
This  figure  establishes  the  attenuation  capability  of  the  ferrite  layer  that  can  be  achieved 
by  tuning  the  permeability.  In  Figure  6  we  graph  the  real  and  imaginary  parts  of  H as 
a  function  of  the  frequency  for  the  incident  angle  0  =  45°.  This  reveals  the  phase  shift 
capability  of  the  controlled  ferrite  layer. 


Figure  5:  Attenuation  Capability 


5  Optimization  of  material  parameters  of  a  coated  air¬ 
foil 

Having  discussed  the  feasibility  of  tuning  dielectric  and  magnetic  properties  of  a  coating  on 
a  perfect  conductor  in  the  previous  sections,  we  turn  in  the  next  several  sections  to  held 
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Figure  6:  Phase  Shift  Property 


calculations  for  a  coated  airfoil  and  demonstrate  our  ability  to  significantly  affect  the  RCS 
by  manipulation  of  the  parameters  e  and  n  in  the  coating. 

5.1  Time-harmonic  Maxwell’s  equation  for  the  transverse  mag¬ 
netic  mode 

We  consider  the  scattering  of  a  perfectly  conducting  airfoil  coated  by  a  layer  of  constant 
thickness.  The  interrogating  electromagnetic  incident  wave  is  assumed  to  be  a  time-harmonic 
and  transverse  magnetic,  more  precisely  TMj.  mode.  Thus,  the  time-harmonic  electric  and 
magnetic  fields  have  the  form 


(5.1) 


where  Ex  is  a  function  of  y  and  z.  We  denote  the  airfoil  by  12  and  the  coating  layer  by  Ox ; 
see  Figure  7. 
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Figure  7:  The  computational  domain,  an  interrogating  wave  E®  with  an  interrogation  angle 
a,  and  an  far  field  pattern  Fa(a  +  n)  of  a  scattered  field  E ^  to  the  direction  a  +  n. 

We  decompose  the  total  electric  field  Ex  to  the  incident  field  E^  and  the  scattered  field 
E^\  that  is,  Ex  =  Ex  ^  +  E^\  Outside  the  coating  layer  the  interrogating  plane  wave  with 
an  interrogation  angle  a  =  (n/2  —  0O)  is  given  by  Ex\y,  z )  =  el^kyy+kzZ\  where  ky  =  kcosa, 
kz  =  k  sin  a,  k  =  2n/\  =  uj/cq  is  the  wave  number,  and  A  is  the  wavelength.  We  do  not 
have  an  explicit  formula  for  the  incident  field  E^’  inside  the  coating  layer  fli.  However  our 
following  formulation  is  constructed  in  such  a  way  that  it  gives  the  same  total  field  Ex  in 
M2  \  0  and  the  same  scattered  field  E ^  in  M2  \  (f^  U  Cl)  regardless  how  the  incident  field 
E is  chosen  in  the  coating  layer  Hi.  Thus,  we  can  use  the  same  expression  el^kyV+kzZ^  for 
the  incident  field  also  in  the  coating  layer. 


By  eliminating  the  magnetic  field  from  Maxwell’s  equation  and  substituting  the  time-harmonic 
electric  field  E  of  the  form  given  in  (5.1)  to  the  resulting  equation,  we  obtain  the  following 
Helmholtz  equation 


V  •  f  -V4r)l  +  ew2^r)  =  -V  ■  ( - V£«  )  -  eu2E «  m2 

F  ) 


Ex  =  Et>  +  Ef  =  0 
1  dEJ  =  [BJ  =  0 


L/i  dn 


in  M2  \  fl 
on  d Cl 

on  <911 1  \  d Cl 


lim  \fr 


r)  ph) 

Hr  ~  **tr)  I  -  «, 


(5.2) 


where  [  •  ]  denotes  the  jump  and  n  is  a  normal  direction  of  the  surface  dCl\  \  dCl.  The  far  field 
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behavior  of  the  scattered  field  described  by  Maxwell’s  equation  satisfies  the  Silver-Miillcr 
radiation  condition  [6].  For  the  time-harmonic  TMX  mode  this  condition  reduces  to  be  the 
Sommerfeld  radiation  condition  given  by  the  limit  in  (5.2).  The  material  permittivity  e  and 
permeability  y,  are  piecewise  constant  functions  defined  by 

,  x  f  ere0,  (y,  z)  e  &i  f  HrHo,  (y,  z )  G 

e{y,z)  =  i  •  and  y{y,z)  =  < 

{  e0,  otherwise  f  /i0 ,  otherwise. 


5.2  Far  field  pattern 

The  far  field  pattern  F  :  [0,  27r]  — >  C  describes  the  intensity  and  phase  of  the  scattered  field 
far  away  from  the  scatterer  ([5],  p.340).  ft  can  be  defined  as 


Fa(0)  =  lim  (Vtoricr  e~i{kr+n/A)  E%\r  cos0,rsin0))  , 


■>oo  \ 


(5.3) 


where  we  have  added  the  subscript  a  to  denote  the  interrogation  angle. 

5.3  Backscatter  reduction 

In  this  optimization  problem,  we  want  to  find  constant  material  parameters  er  and  /ir  so 
that  the  intensity  of  the  backscattered  wave  is  minimized  over  a  given  sector  [ao,ai\.  The 
objective  function  is  the  integral 


J(er,fj,r)=  /  \Fa(ot  +  7r)|2  da. 


(5.4) 


'  OL  0 


The  minimization  problem  is  given  by 


min  J(er,yr), 

(er,/ir)sQ 

where  Q  is  the  set  of  admissible  material  parameters. 


(5.5) 


This  objective  function  corresponds  to  a  situation  where  the  same  radar  is  illuminating  and 
detecting  the  scattered  wave.  The  interrogation  angle  of  the  interrogating  wave  varies  within 
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the  interval  [ao,«i]  and  the  formulation  (5.4)  corresponds  to  an  assumption  of  a  uniform 
distribution  on  possible  angles  of  interrogation. 

We  use  the  NAG  Fortran  library’s  [17]  E04UCF  implementation  of  a  sequential  quadratic 
programming  (SQP)  method.  This  is  a  gradient  based  optimization  method  which  approxi¬ 
mates  the  gradient  using  finite  differences. 


6  Approximation 

6.1  Truncation  of  domain  and  variational  formulation 


For  the  discretization  of  (5.2),  we  restrict  the  problem  to  a  rectangular  domain  ff  and 
impose  a  second-order  absorbing  boundary  condition  [1]  on  the  artificial  boundary  dll  to 
approximate  the  the  Sommerfcld  radiation  condition.  Now  the  scattered  field  E 4  satisfies 
the  following  equations 


V  •  (  -v4r)  I  +  ew24r)  =  -V  •  (  -vM<} 


euj2E^ 


1  dEr 


F  dn 


i  d2E{£) 


r)  F ^ 

J  x _ ikE (r) _ 

dn  x  2k  ds2 


4r)  = 

=  [Ex\  —  0 
=  0 


dE. 


(r) 


ds 


ik-Eir)  =  0 
2 


in  fl  \  Cl 
on  dfl 

on  <9Qi  \  dfl  (g_x) 
on  tin 
at  C, 


where  n  and  s  denote  the  normal  and  tangential  directions  of  the  boundary  <911,  respectively, 
and  C  is  the  set  of  the  corner  points  of  ff. 


The  variational  formulation  of  (6.1)  is  : 
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6.2  Finite  element  approximation 

The  finite  element  approximation  is  implemented  using  linear  elements.  The  mesh  is  con¬ 
structed  from  two  uniform  triangular  meshes.  The  finer  mesh  is  for  the  coating  layer  and, 
the  coarser  is  for  the  exterior  domain  outside  the  coating.  The  mesh  step  sizes  are  chosen  in 
such  a  way  that  the  number  of  nodes  per  wavelength  is  approximately  the  same  in  the  air 
and  coating.  The  finer  mesh  is  locally  fitted  to  the  surfaces  of  the  obstacle  and  the  coating 
layer.  The  local  fitting  is  done  using  the  algorithm  [4]  with  slight  modifications.  Between 
the  meshes,  there  is  a  layer,  which  fits  the  meshes  together  in  a  conforming  way.  An  example 
of  a  part  of  mesh  is  shown  in  Figure  8. 


Figure  8:  A  magnified  view  of  the  mesh  for  a  coated  NACA0012  airfoil. 


After  the  discretization  of  the  variational  formulation  (6.2),  we  obtain  a  system  of  linear 
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equations 


Ax  =  b,  (6.3) 

where  the  matrix  A  is  a  symmetric  non-Hermitian  complex  matrix.  The  vector  x  contains 
the  nodal  values  of  the  scattered  held  E ^  and  the  vector  b  corresponds  the  right-hand  terms 
in  (6.2).  The  resulting  systems  of  linear  equations  are  solved  using  an  iterative  method  which 
combines  fictitious  domain  and  domain  decomposition  methods  [8]. 


The  far  held  pattern  Fa(a  +  tt)  in  (5.3)  is  computed  as  a  surface  integral  [6]  of  the  computed 
near  held  E(xr'1  and  its  hux.  Our  particular  implementation  of  the  computations  is  described 
in  [7].  We  need  to  evaluate  the  far  held  pattern  Fa(a  +  tt)  in  (5.3)  using  the  computed  near 
held  E^\  Our  particular  implementation  of  the  computations  is  described  in  [7].  The  basic 
idea  of  this  procedure  is  the  following.  Let  E ^  be  the  harmonic  extension  of  E ^  from 
fl\  (hUhj)  to  M2.  We  obtain  this  extension  as  a  byproduct  of  our  solution  procedure.  Then 
at  the  point  rj  outside  0  U  Oi  the  scattered  held  E ^  is  given  by 


£«(,,)  =  f  (A £$■>«)  +  k2EM(0)  i(r,,£)d£,  (6.4) 

where  $(77,  £)  is  the  fundamental  solution  of  the  homogeneous  Helmholtz  equation  in  M2 
given  by  the  Hankel  function  $(77,  £)  =  -H^\k\r]  —  £|).  By  using  Green’s  formula  one  can 
show  that  (6.4)  is  equivalent  to 


Eir\v)  = 


Ei:Ha)aiM 


dE{l\ 


cr 


1  an 


dn 


dn 


$(77,  cr)  )  da 


which  is  a  more  traditional  expression  for  E)p  (rj)  [6].  The  far  held  pattern  is  obtained  by 
hrst  discretizing  the  Helmholtz  operator  and  integral  in  (6.4)  and  then  taking  the  limit 


Fa(a  +  7 r)  =  lim  ( V 8nkr  e  *(fcr+,r/4)  E^\r  cos(ck  +  tt),  r  sin(a  +  tt)) 


At  the  discrete  level  the  previous  procedure  reduces  to  the  evaluation  of  a  sum  of  exponential 
functions.  This  can  be  performed  easily  and  fast. 
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We  have  compared  numerical  results  computed  using  the  proposed  method  and  the  same 
implementation  with  test  cases  presented  on  a  workshop  in  Oxford  [18].  The  two-dimensional 
test  cases  did  not  include  a  coated  NACA  airfoil,  but  they  did  include  a  similar  coated  ogive 
(the  intersection  of  two  non-concentric  disks).  We  computed  numerical  results  for  10  test 
cases.  All  our  numerical  results  and  especially  the  radar  cross  sections  were  in  very  good 
agreement  with  the  majority  of  the  results  presented  in  the  workshop  for  each  test  case. 
Because  of  this,  we  expect  our  results  to  be  accurate  also  for  the  coated  NACA  airfoils. 

7  Numerical  experiments 

In  our  experiments,  we  minimize  the  backscatter  by  a  coated  NACA0012  airfoil.  The  length 
of  the  airfoil  is  one  unit  without  coating  and  the  trailing  edge  of  the  perfectly  conducting 
material  is  at  the  origin.  We  minimize  the  backscatter  for  the  interrogation  angles  in  the 
sector  [ao^i]  =  [0, 7r] .  We  considered  two  wavelengths,  A  =  1/4  and  A  =  1/10.  The 
thickness  of  the  coating  is  A/ 10. 

We  consider  first  the  lower  frequency  experiments  with  a  four  wavelength  long  airfoil.  The 
computational  domain  is  [—1.5, 0.5]  x  [—0.6, 0.6].  Our  discretization  has  20  nodes  per  wave¬ 
length  in  the  ambient  medium  leading  to  a  triangulation  with  22157  nodes  and  43158  ele¬ 
ments.  A  magnified  view  of  the  mesh  at  the  trailing  edge  is  shown  in  Figure  8.  We  perform 
several  optimizations  with  different  box  constraints  for  the  real  and  imaginary  parts  of  er  and 
fir.  The  results  of  these  optimizations  are  given  in  Table  1.  A  surface  plot  of  the  objective 
function  is  shown  in  Figure  9.  The  radar  cross  sections  defined  by 

Fa(a  +  vr)|-j 

are  shown  for  two  optimized  materials  in  Figure  10.  Corresponding  reflected  held  inten¬ 
sities  for  comparison  between  the  no  coating  layer  case  and  the  optimized  complex-valued 
parameters  case  of  Figure  10  are  depicted  in  Figures  11  and  12  for  an  angle  of  interrogation 


RCS(a)  =  10  log10  (  — 
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1 

0 

0 

0 

324.340 

€r 

1 

3.92 

10 

0 

0 

0 

1 

1 

10 

0 

0 

0 

318.738 

6r 

1 

6.52 

10 

0 

3.14 

10 

fir 

1 

1 

1 

0 

0 

0 

2.918 

6r 

1 

7.75 

10 

0 

0.80 

10 

fir 

1 

5.41 

10 

0 

2.51 

10 

0.083 

Table  1:  The  results  of  material  parameter  optimizations  for  the  case  A  =  1/4. 


Real  part 
min  opt  max 

Imaginary  part 
min  opt  max 

J (6r ,  /^r) 

6r 

1 

1 

1 

0 

0 

0 

hr 

1 

1 

1 

0 

0 

0 

806.036 

6r 

1 

3.84 

10 

0 

0 

0 

hr 

1 

1 

10 

0 

0 

0 

802.695 

6r 

1 

6.58 

10 

0 

3.14 

10 

fir 

1 

1 

1 

0 

0 

0 

2.901 

6r 

1 

4.54 

10 

0 

2.99 

10 

hr 

1 

4.99 

10 

0 

2.92 

10 

0.234 

Table  2:  The  results  of  material  parameter  optimizations  for  the  case  A  =  1/10. 


a  =  7t/4. 


The  computational  domain  for  the  higher  frequency  experiments  with  a  ten  wavelengths 
long  airfoil  is  [—1.3,  0.3]  x  [—0.4,  0.4],  Again  our  discretization  has  20  nodes  per  wavelength 
in  the  ambient  medium  leading  to  a  triangulation  with  65551  nodes  and  128530  elements. 
The  results  of  optimizations  with  different  box  constraints  are  given  in  Table  2.  Radar  cross 
sections  for  two  optimized  materials  are  shown  in  Figure  13.  Corresponding  reflected  field 
intensities  for  comparison  between  the  no  coating  layer  case  and  the  optimized  complex¬ 
valued  parameters  case  of  Figure  13  are  depicted  in  Figures  14  and  15  for  an  angle  of 
interrogation  a  =  7t/4. 
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0  10 

Figure  9:  The  objective  function  J(er,/ir)  for  complex-valued  er  and  fir  —  1  for  the  case 
A  =  1/4. 


Figure  10:  The  RCS  for  optimized  complex-valued  material  parameters  (er  =  7.75  +  0.807; 
Hr  =  5.41  +  2.517;  solid  line),  optimized  complex- valued  permittivity  (er  =  6.52  +  3. 147; 
Hr  —  1;  dotted  line),  and  for  no  coating  (er  =  1;  fir  =  1;  dashed  line)  for  the  case  A  =  1/4. 
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Figure  11:  The  reflected  field  E^fr)  for  no  coating  (er  =  1;  /ir  =  1)  for  the  case  A  =  1/4  and 
angle  of  interrogation  a  =  7t/4. 


Figure  12:  The  reflected  field  E^'1  for  optimized  complex-valued  material  parameters  (er  = 
7.75  +  0.807;  /ir  =  5.41  +  2.517)  for  the  case  A  =  1/4  and  angle  of  interrogation  a  =  7r/4. 
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Figure  13:  The  RCS  for  optimized  complex-valued  material  parameters  (er  =  4.54  +  2.997; 
/ir  =  4.99  +  2.92;  solid  line),  optimized  complex- valued  permittivity  (er  =  6.58  +  3. 147; 
Hr  =  1;  dotted  line),  and  for  no  coating  (er  —  1 ;  jur  —  1;  dashed  line)  for  the  case  A  =  1/10. 


26 


Figure  14:  The  reflected  field  E ^  for  no  coating  (er  =  1;  /ir  =  1)  for  the  case  A  =  1/10  and 
angle  of  interrogation  a  =  7t/4. 


Figure  15:  The  reflected  field  E^'1  for  optimized  complex-valued  material  parameters  (er  = 
4.54  +  2.99i;  /jr  =  4.99  +  2.92f)  for  the  case  A  =  1/10  and  angle  of  interrogation  a  =  it /A. 
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8  Summary  and  conclusions 


In  summary,  the  efforts  reported  on  in  this  paper  are  an  important  first  step  in  developing  an 
attenuation  or  anti- interrogation  technology.  We  considered  first  the  question  of  the  feasibil¬ 
ity  of  reduction  in  reflected  electromagnetic  waves  from  a  planar  layered  coating  on  a  perfect 
conductor.  We  demonstrated  that  even  under  uncertainty  of  the  interrogating  wavelengths 
(frequencies),  one  can  achieve  reduction  of  the  reflection  coefficient  through  optimizing  the 
dielectric  permittivity  in  a  coating  layer.  We  then  considered  a  Maxwell  equation  based  for¬ 
mulation  for  a  composite  ferromagnetic-ferroelectric  device  built  and  experimentally  tested 
by  How  and  Vittoria.  We  derived  the  pertinent  reflection  field  equations  for  time  harmonic 
interrogating  TMX  mode  plane  waves  and  showed  that  substantial  control  of  reflected  waves 
(both  magnitude  and  phase)  can  be  obtained  by  tuning  the  magnetic  permeability  of  a  ferrite 
layer.  We  next  turned  to  a  non-planar  geometry,  in  this  case  a  2-D  airfoil,  with  a  coating 
layer  wherein  both  the  dielectric  permittivity  e  and  the  magnetic  permeability  //  can  be 
optimized.  Allowing  a  uniform  uncertainty  on  the  interrogating  signal  angles,  we  use  com¬ 
putational  methods  from  2-D  scattering  theory  (the  Helmholtz  equation  with  Sommcrfcld 
far  field  radiation  conditions)  to  verify  that  significant  reduction  in  the  far  field  reflection 
can  be  obtained  by  optimal  choice  of  e  and  /i. 

All  of  the  investigations  discussed  in  this  paper  were  pursued  under  an  active  design  scenario 
and  did  not  allow  for  online  adaptivity  of  the  coating  layers.  Future  investigations  of  great 
interest  include  the  feasibility  of  combining  the  formulations  in  this  paper  with  real  time 
sensing  and  adaptive  (feedback)  control  of  coatings  such  as  those  described  above  to  develop 
an  active  control  attenuation  capability. 
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